RANDOM PERMUTATION MATRICES UNDER THE 
GENERALIZED EWENS MEASURE 

CHRISTOPHER HUGHES, JOSEPH NAJNUDEL, ASHKAN NIKEGHBALI, 

AND DIRK ZEINDLER 



Abstract. We consider a generalisation of the Ewens measure for the 
symmetric group, calculating moments of the characteristic polynomial 
and similar multiplicative statistics. In addition we study the asymp- 
totic behaviour of linear statistics (such as the trace of a permutation 
matrix or of a wreath product) under this new measure. 



Contents 

1. Introduction ^ 
Acknowledgements ^ 

2. The generalised Ewens measure, generating series and singularity 
analysis ^ 

2.1. Some combinatorial lemmas and the generalised Ewens measure ^ 

2.2. Generating functions and singularity analysis ^ 

3. Moments of multiplicative class functions 13 

3.1. Multiplicative class functions 13 

3.2. Generating functions for 15 

3.3. Asymptotic behaviour of the moments 16 

4. Traces |18 

4.1. Traces of permutation matrices 19 

4.2. Traces of functions 21 

5. Wreath product, traces and the generalised Ewens measure 23 

6. Diverging variance for the classical Ewens measure 31 
References 37 



Date: September 26, 2011. 



1 



2 



C. HUGHES, J. NAJNUDEL, A. NIKEGHBALI, AND D. ZEINDLER 



1. Introduction 

The trace of a matrix is one of the most natural additive class functions 
associated to the spectra of a matrix. Traces of unitary matrices chosen ran- 
domly with Haar measure have been much studied, for example by Diaconis- 
Shahshahani |^], Diaconis-Evans and Rains p^] using methods from rep- 
resentation theory; by Diaconis and Gamburd M using combinatorics; and 



using methods from mathematical physics by Haake et al. 

Another natural class function, this time multiplicative, is the characteristic 
polynomial, and the distribution of characteristic polynomials of random 
unitary matrices has been studied by many authors, including Keating- 
Snaith [17] and Hughes, Keating and O'Connell 



From the characteristic polynomial one can find the number of eigenval- 
ues lying in a certain arc (since the underlying matrix is unitary, all the 
eigenvalues lie on the unit circle). The problem of studying the number of 
eigenvalues lying in an arc was studied by Rains pi] and Wieand ]26], who 



found a very interesting correlation structure when multiple arcs were con- 
sidered, and Hughes, Keating and O'Connell ]jl^ who made the connection 
with characteristic polynomials. 

One of the reasons for such an extensive study into random unitary matrices 
and their spectra is that the statistical distribution of the eigenvalues is 
expected to have the same behaviour as the zeros of the Riemann zeta 



function (see Montgomery [21] and Keating-Snaith [|l] 

The statistics of the distribution of spectra of infinite subgroups of the uni- 
tary group, such as the symplectic and orthogonal groups are also expected 
to model families of other L- functions (Keating-Snaith [^]), and the distri- 
bution of traces for these groups have been studied, frequently in the same 
papers. 

However, the statistics of the spectra of finite subgroups of the unitary 
group, such as the permutation group, is not so well studied, though there 
are many results known. 

Wieand ]25] studied the number of eigenvalues of a uniformly random per- 
mutation matrix lying in a fixed arc, and Hambly, Keevash, O'Connell and 



Stark ]12] found corresponding results for the characteristic polynomial, 
making the same connection between the characteristic polynomial and the 
counting function of eigenvalues. 

In all these cases, the permutation matrices were chosen with uniform mea- 
sure and the results were similar to those found for the full unitary group 
with Haar measure. However, there were some significant differences primar- 
ily stemming from the fact that the full unitary group is rotation invariant, so 
the characteristic polynomial is isotropic. The group of permutation matrix 
is clearly not rotation invariant, and the distribution of the characteristic 
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polynomial depends weakly on the angle of the parameter. Most results 
require the angle to have the form Ittt with r irrational and of finite type. 
(It is worth pointing out that those angles which are not of this type, have 
Hausdorff dimension zero). More recently Ben Arous and Dang (|Q) have 
extended some of the results of Wieand to more general measures together 
with new observations very specific to permutation matrices. In particular 
they prove that the fluctuations of smooth linear statistics (with bounded 
variance) of random permutation matrices sampled under the Ewens mea- 
sure are asymptotically non-Gaussian but infinitely divisible. 

A permutation matrix has the advantage that the eigenvalues are determined 
by the cycle-type A = (Ai, • • • , A^) of the corresponding permutation. This 
allows to write down in many situations explicit expression for the studied 
object. For example, the characteristic polynomial of an x permutation 
matrix M is 

ZNix) = det{I-xM) = YI{1- xe^°^) = ~ ^^'") = 11 " ^'') 

j=l m=l k=l 

(1.1) 

where e*"^ , • • • , e*°^ are the eigenvalues of M and Ck is the number of cycles 
of length k (i.e. the number of j such that Xj = k). 

Equation ( |1 . ID has been used by Dehaye and Zeindler to introduce mul- 
tiplicative class functions (i.e. invariant under conjugation) associated to 
function /. These functions have the form 

H^^(/)(^)=n/M ' (1-2) 

k=l 

and generalise the characteristic polynomial, which is the case f{x) = (1— x). 
This nice generalisation is not possible for the unitary group because the 
eigenvalues lack the structure of eigenvalues of permutation matrices that 
allows ( |1 . 1[ ) to hold. The most natural analogue for unitary matrices is 
Heine's identity which connects Haar averages of multiplicative class func- 
tions with determinants of Toeplitz matrices. 

Equation (^]^) can also be used to express linear statistics or traces of func- 
tions of random permutation matrices in terms of the cycle counts. More 
precisely, if we identify a random permutation matrix M with the permuta- 
tion a it represents, we have: 

Definition 1.1. Let F : ^ C be given. We then define the trace of F to 
be the function Tr(i<') : ©at — t- C with 

N 

Tr{F)ia):=Y,nook) (1-3) 

k=l 

where {oJk)^^i are the eigenvalues of a with multiplicity. 
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Observe that when F{x) = x'^, we have TT{F){a) = Tr(cr°'), and this justifies 
the use of the terminology trace. The trace of a function is also referred to 
as a linear statistic on ©at. 

Now it follows immediately from equation ( |1.1D that: 

Lemma 1.2. Let F : ^ C and a £ &iy with cycle type X be given. We 
then have 

N 

TiiF){a) = Y,kCkAkiF) (1.4) 

k=l 

With Ak(F) :=iELi^(e'"*'"/')- 

The reason why the expression in ( |l.lD and (|l.2D are useful is that many 
things are known about the cycle counts Ck- For example, the cycle counts 
Ck converge weakly to independent Poisson random variables Pk, with mean 
1/k. Also useful in this context is the Feller-coupling. This allows in many 
situations to replace by P^- Several details on the cycle counts and the 
Feller coupling can be found in the book by Arratia, Barbour and Tavare 

These results all concern the uniform measure, where each permutation has 
weight 1/A^!, or the Ewens measure, where the probability is proportional to 
the total number of cycles, 9^^"^^ / {Nlhjsi), with /i^v the required normalisation 
constant (the case 9=1 corresponding to the uniform measure) . A common 
ingredient in all the above cited works is the use of the Feller coupling (details 
can again be found in the book by Arratia, Barbour and Tavare Q) and 
some improvements on the known bounds for the approximation given by 
this coupling (see Q, Section 4). 

In recent years, there have been many works in random matrix theory aim- 
ing at understanding how much the spectral properties of random matrices 
depend on the probability distributions of its entries. Here a similar ques- 
tion translates into how are the linear and multiplicative (i.e. multiplicative 
class functions) statistics affected if one considers more general probabil- 
ity distributions than the Ewens measure on the symmetric group? The 
Ewens measure can be naturally generalised to a weighted probability mea- 
sure which assigns to the permutation matrix M (i.e. to the associated 
permutation) the weight 

1 ^ C 

where /i at is a normalisation constant. The Ewens measure corresponds to 
the special case where 6k = 9 is a constant. This measure has recently ap- 
peared in mathematical physics models (see e.g. Q and Q) and one has 
only recently started to gain insight into the cycle structures of such random 
permutations. One major obstacle with such measures is that there exists 
nothing such as the Feller coupling and therefore the classical probabilistic 
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arguments do not apply here. In a recent work, Nikeghbali and Zeindler 
[p3| propose a new approach based on combinatorial arguments and singu- 
larity analysis of generating functions to obtain asymptotic expansions for 
the characteristic functions of the C^s as well as the total number of cycles, 
thus extending the classical limit theorems (and some distributional approx- 
imations) for the cycle structures of random permutations under the Ewens 
measure. In this paper we shall use the methods introduced in namely 
some combinatorial lemmas, generating series and singularity analysis, to 
study linear and multiplicative statistics for random permutation matrices 
under the general weighted probability measure. In fact we shall consider 
the more general random matrix model obtained from the wreath product 
I &M (see for instance [^]): this amounts to replacing the I's in the 
permutation matrices by independent random variables taking values in the 
unit circle S^. It should be noted that many groups closely related to Sat 
exhibit such matrices, for instance the Weyl group of S0{2N). 

More precisely this paper is organised as follows: 

In Section Q we fix some notation and terminology, recall some useful com- 
binatorial lemmas together with some results of Hwang (and some slight 
extensions) on singularity analysis of generating functions. In particular we 
shall introduce two relevant classes of generating functions according to their 
behaviour near singularities on the circle of convergence. In this article, we 
shall state our theorems for random matrices under the generalised Ewens 
measures for which the generating series of {9k)k>i is in one of these two 
classes. 

In Section |^ we study the multiplicative class functions associated to a func- 
tion / and obtain the asymptotic behaviour of the joint moments. In par- 
ticular we extend earlier results of [^], [28] and of [12| on the characteristic 



polynomial of uniformly chosen random permutation matrices. 

In Section ^ we focus both on the traces of powers and powers of traces 
which are classical statistics in random matrix theory. In fact we prove 
more generally that the fluctuations of the linear statistics for Laurent poly- 
nomials are asymptotically infinitely divisible (they converge in law to an 
infinite weighted sum of independent Poisson variables). We also establish 
the convergence of the integer moments of linear statistics for functions of 
bounded variation together with the rate of convergence. 

In Section ^ we consider the more general model consisting of the wreath 
product IGn and study the linear statistics for general functions F in 
(1.3). In such models, the I's in the permutation matrix are replaced with 



(zj)i<j<jv which are i.i.d random variables taking their values on the unit 



circle S^. In this framework. Lemma 1.2 can be naturally extended (see 
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Lemma 5^) and the quantity 

Afc(F,z):=i ^F(a;). (1.5) 

naturally appears in our technical conditions. Under some conditions on rate 
of convergence to of the L^-norm of Afc(F, z), and some assumptions on 
the singularities of the generating series of {6k)k>i, we are able to compute 
the asymptotics of the characteristic function of Tr{F) with a good error 
term. From these asymptotics we are able to compute the fluctuations of 
Tt{F). We also translate our conditions in terms of the Fourier coefficients 
of F, where F has to be in some Sobolev space H^. 

In Section |6| we still work within the framework of the wreath product S^IGn 
and consider the case where the variance of Tr(F) is diverging. This time 
we restrict ourselves to the Ewens measure since our methods do not seem 
to apply in this situation. We hence go back to probabilistic arguments (i.e. 
use the Feller coupling) to prove that under some technical conditions on 
F, the fluctuations are Gaussian. In fact we essentially adapt the proof by 
Ben Arous and Dang ([^) to this more general situations. Nonetheless our 
theorem is also slightly more general in that it applies to a larger class of 
functions F. 
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2. The generalised Ewens measure, generating series and 

singularity analysis 

In this section we fix notation and we recall some facts about the symmetric 
group and generating functions, as well as the main results from singular 
analysis (we also provide some variants and extensions for the purpose of 
this paper). Our presentation follows closely |2^ ]. 



2.1. Some combinatorial lemmas and the generalised Ewens mea- 
sure. We present in this section some basic facts about &n and then define 
the generalised Ewens measure. We give here only a very short overview 
and refer to |l| and [20| for more details. 



2.1.1. Conjugation classes and functions on &n- We first take a closer look 
at the conjugation classes of the symmetric group Sat (the group of all 
permutations of a set of objects). We only need to consider the conjuga- 
tion classes since all probability measures and functions considered in this 
paper are invariant under conjugation (i.e. they are class functions). It is 



RANDOM PERMUTATION MATRICES UNDER THE GENERALIZED EWENS MEASURE 



well-known that the conjugation classes of &N can be parameterised with 
partitions of N. 

Definition 2.1. A partition X is a sequence of non-negative integers Ai > 
A2 > • • • eventually trailing to O's, usually omitted. The size of the partition 
is |A| := Ylim^rn- We call A a partition of N if |A| = A^, and this will he 
denoted by \ \- N. The length of A is the largest I such that A^ / 0. 

Let (T € ©AT be arbitrary. We can write a = g\ - ■ ■ oi with 0"^, 1 < m < i, 
disjoint cycles of length A^. Since disjoint cycles commute, we can assume 
that Ai > A2 > • • • > A^. We call the partition A = (Ai,A2,-- - ,Xi) the 
cycle-type of a. We write Cx for the set of all a G ©at with cycle type A. 
One can now show that two elements a,T £ ©tv are conjugate if and only if 
a and r have the same cycle- type and that Cx are the conjugation classes of 



©AT. Since this is well-known, we omit the proof and refer to |20| for more 
details. 

Definition 2.2. Let o" G ©at be given with cycle-type A. The cycle numbers 
Ck are defined as 

Ck = Ck{a):=#{m:\ra = k} (2.1) 
and the total number of cycles T{a) is 

N 

T{a):=^Ck. (2.2) 
fc=i 

The functions Cjt(cr) and T{a) depend only on the cycle type of a are thus 
class functions. Clearly T(cj) equals ^(A), the length of the partition corre- 
sponding to a. 

All expectations in this paper have the form X^o-66]v ^('^) ^'-'^ ^ class 
function u. Since u is constant on conjugation classes, it is more natural 
to sum over all conjugation classes. We thus need to know the size of each 
conjugation class. 

Lemma 2.3. We have 

I© I ^ 
|Ca| = ^—^ with zx := JJ '^^"^kl (2.3) 

k=i 

with Ck defined in ( |2.1| ), and 

^ E = E ^-(Ca) (2.4) 

' creejv AhAf ^ 

for a class function u : <3n ^ C 



Proof. The first part can be found in [20| or in chapter 39]. The second 
part follow immediately from the first part. □ 
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2.1.2. Definition of the generalised Ewens measures. We now define the gen- 
eralised Ewens measures. 

Definition 2.4. Let Q = {9k)^=i be a sequence of strictly positive numbers. 
We define for a £ &j\f with cycle-type A 

771=1 k = l 

with Hn = h]\j{Q) a normalisation constant and Hq := 1. 

The second equahty in (|2.5|) fohows immediately from the definition of Ck 



(Definition 2.2). The uniform measure and the Ewens measure are special 
cases, with 9i^ = \ and = 9 a, constant respectively. 

We now introduce two generating functions closely related to Pe- 



oo /, / oo 



■= E f Ge(t) := exp E f (2.6) 

k=l \k=l / 

At the moment, (70 (t) and GQ{t) are just formal power series, however we 
will see in Section |2.2| that 



Ge{t) = h^t^. (2.7) 

Af=0 

where the h n are given in Definition 12.41. 



2.2. Generating functions and singularity analysis. The idea of gen- 
erating functions is to encode information of a sequence into a formal power 
series. 

Definition 2.5. Let {gN)j>^^^ be a sequence of complex numbers and define 
the (ordinary) generating function of the sequence as the formal power series 



G{t) = 9Nt''. (2.. 
We define [t^] [G] to he the coefficient of t^ in G{t), i.e. [t^] [G] := qn, 



The reason why generating functions are useful is that it is often possible to 
compute the generating function without knowing gi\j explicitly. 

The main tool in this paper to calculate generating functions is: 
Lemma 2.6. Let {am)meN be a sequence of complex numbers. Then 

with the same z\ as in Lemma W^. 
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// any one of the sums in is absolutely convergent then so are the 

others. 

Proof. The first equality follows immediately from the definition of Ck • The 
proof of the second equality in (|2.9| ) can be found in or can be directly 
verified using the definitions of zx and the exponential function. The last 
statement follows with dominated convergence. □ 

We now use this lemma to prove the identity given in ( |2.7D . The constant 
h]y in ( |2.5| ) is chosen so that P0[o"] is a probability measure on 6n- It thus 
follows that 

■ o-e6jv k=l \\-N ^ m=l 

It now follows with Lemma p.6| that 



E ^^^'^ = E f n = exp 5: f = Geit), (2.11) 

Af=0 A m=l \k=l J 



which proves (^ 

Corollary 2.7. In the special case of the Ewens measure, when Oj. is a 
constant 9, say, we have 



Ge{t)=Y.hMt'' = {l-t)-'. (2.12) 



Af=0 



From this it immediately follows that = (— 1) 



Given a generating function G{t), a natural question is: what is the co- 
efficient of t^ and what is the asymptotic behaviour of this coefficient as 
N — )• cxD? If G{t) is holomorphic near then one can use Cauchy's integral 
formula to do this. But it turns out that it is often difficult to compute the 
integral exactly, but we will now see that one can nevertheless extract the 
asymptotic behaviour of the coefficient when G{t) has a special form. 

Definition 2.8. Given R> r and < < |, let 

Ao = Ao(r, iJ, (^) = {z G C : |z| < i?, z / r, I arg(z - r)| > 0} (2.13) 



Definition 2.9. Let r > 0, ■!? > and a complex constant K he given. We 
say that a function g(t) is in J-{r, "d, K) if there exists R > r and < (f) < ^ 
such that g{t) is holomorphic in Ao{r,R,(p), and 

g{t) = mog (^-^-^jj-^ + K + O {t - r) (2.14) 

as t ^ r with t G Aq (r, i?, (/>) . 
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The following theorem, proven by Hwang in |[T^, gives the asymptotic be- 
haviour of the coefficient of t'^ for certain special generating functions 

Theorem 2.10 (Hwang Q). Let g{t) £ T{r,'&,K), and let 5(t) be holo- 
morphic in \t\ < r. Set G{t,w) = e^^^^^ S{t), then 

[,«„G,,„),.£!^!^f^,0fll) ,2.15) 



uniformly for bounded complex w. 

Remark. The idea of the proof is to take a suitable Hankel contour and 
to estimate the integral over each piece. The details can be found in |14, 
chapter 5]. 

Remark. One can compute lower order error terms if one has more terms in 
the expansion of g{t) near r. 

As a first simple application of this result, we compute the asymptotic be- 
haviour of h]\r for the generalised Ewens measure if gQ{t), as defined in (2.6), 
is in J^{r, "i?, K). 

Lemma 2.11. Let ge{t) S T{r,'i!}, K). We then have 



e 



Proof. We have proven in ( |2.11| ) that Yl'N=o^Nt^ = exp(ge(t)). We thus 
can apply Theorem 2.1C with g{t) = geit), w = 1 and S{t) = 1. □ 

Remark. For the Ewens measure, when Q is the constant sequence {0)^i, 
we have ge{t) G T{1,6,0) and thus = ^Ylef (l + ^ {jf))- However, in 
this special case, one can do much more since Hn is known to equal {^~^^~^) ■ 
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Essentially, one can think of Hwang's result as concerning functions with a 



solitary singularity at t = r. In Section 3.2 we will need a version of this 
theorem with multiple singularities that we now state. 

Definition 2.12. Let ^ = (Cj)f=i with / for i 7^ j, and with = r 
(so the are distinct points lying on the circle of radius r). Let R > r and 
let < 4> < ^, then set 

d 

Arf(r,i?,(/.,e) := P|{z G C : |z| < i?, z ^ ^i, |arg(z - Ci) - arg(ei)| > 0} 
1=1 

(2.17) 




Figure 2. Illustration of the domain Ad{r, R, cp,^) 

Definition 2.13. Let ^ = (i?j)^^^ and K = (Ki)'^^-^ be two sequences of 
complex numbers, and let r > 0. We say a function g{t) is in J^{r,'d,K) 
if there exists R > r and < (f) < ^ such that g{t) is holomorphic in 
Arf(r, R, (p,C), and for each i = 1, . . . , d, 

g{t) = 'dilog (^-^-^^ + Ki + O {t - ^i) (2.18) 

as t ^ with t G Ad{r, R, (p, ^) . 



Theorem p.lOj generalises to 

Theorem 2.14. Let g G T{r,'&,K), and let S{t) be holomorphic in t for 
\t\ < r. Set G{t,w) = e'"9ii)s{t). We have 

n,c,,.),.|:^(^.o(i)) (...) 
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uniformly for bounded w. 



Sketch of the proof. The proof is combination of the proof of a multiple sin- 
gularities theorem in [10, Section VI. 5] and the proof of Theorem 2.10| . More 
precisely, we apply Cauchy's integral formula with the curve C illustrated 
in Figure ^, where the radius R of the great circle is chosen fix with R> r, 
while the radii of the small circles are 1/n. 




Figure 3. The curve C 



A straightforward computation then shows that the integral over this curve 
gives ( p. 191) and that the error terms are uniform for bounded w. □ 

In practice the computation of the asymptotic behaviour near the singularity 
is often very difficult, and it is not easy to prove whether a function g[t) is in 
J-{r,'d,K) or not. An alternate approach is to combine singularity analysis 
with more elementary methods. The idea is to write G = G1G2 in a way that 
we can apply singularity analysis on Gi and can estimate the growth rate of 
[t^] [G2]. One then can compute the coefficient [t^] [G] directly and apply 
elementary analysis on it. This method is called the convolution method. 

Definition 2.15. Let i?>0, r>0, 0<7< 1 be given. We say g{t) is in 
eJ^(r, t?,7) if g{t) is holomorphic in \t\ <r with 

g{t)=i)log(r^-^)+goit) (2.20) 
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and 

[t""] N = 0(r-^iV-i-^) (2.21) 

as N ^ oo. 



Theorem 2.16 (Hwang ||T|]). Let g{t) G eJ"(r,i9,7), and let S{t) he holo- 
morphic in \t\ < r. Set G{t,w) = e^^^^^S{t), then 



Afi r^/. . .M _ e^'SoMiV""'-! S{r) 



with 



n = pv f(i^ (2.22) 



O ZTJ Re('u;) > 



Rn{w) = { 1 ,x ^ (2.23) 



uniformly for bounded w. 



This theorem is more general than Theorem 2.10, but the error terms are 
worse. As in lemma 2.11, we can compute the asymptotic behaviour of /i^v 



in the case of the generalised Ewens measures when gQ,{t) G eJ-"(r, -i?, 7). 
Lemma 2.17. Assume that ge{t) G eJ^(r, i?, 7). We then have 

3. Moments of multiplicative class functions 



We extend in this section the results of [|T2[, and ||5[ to the generalised 
Ewens measure Fq. More precisely, we compute the asymptotic behaviour 
of the moments of the characteristic polynomial Zj^{x) and of multiplicative 
class functions [P') with respect to Pe using the methods of generating 
functions and singularity analysis introduced in the previous section. 

3.1. Multiplicative class functions. It is well known that Sat can be 
identified with the group of permutation matrices via 

It is easy to see that this map is an injective group homomorphism. We 
thus do not distinguish between S tv and the group of permutation matrices 
and use for both the notation Sat. It will be always clear from the context 
if it is necessary to consider cr G Sat as a matrix. 

Definition 3.1. Let x G C and a G Sat. The characteristic polynomial of 
a is 

Zn{x) = Zj^{x){a) := det{lN — xa). (3-2) 

It is a standard fact that the characteristic polynomial can be written in 
terms of the cycle type of a. 
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Lemma 3.2. Let a £ 6n be given with cycle type X; then 

ZNix)= Ylil-x^"^), (3.3) 

m=l 

with i{X) the length of the partition A, which is the same as the number of 
cycles T{(t). 

Proof. Since any permutation matrix is conjugate to a block matrix with 
each block corresponding to one of the cycles, and the characteristic polyno- 
mial factors over the blocks, it is sufficient to prove this result in the simple 
case of a one-cycle permutation, where it follows from a simple calculation. 
More explicit details can be found for instance in p^, Chapter 2.2]. □ 



Equation (3^) shows that the spectrum of permutation matrix is uniquely 
determined by the cycle type. We use this as motivation to define multi- 
plicative class functions on ©at. 

Definition 3.3. Let P{x) he a polynomial in x. We then define the multi- 
plicative class function associated to the polynomial P as 

W^iP){x) = W^iP){x){a) := n ^ (^^'") ' (3-^) 

m=l 

For brevity, we simply call this a multiplicative class function. 

It follows immediately that the characteristic polynomial is the multiplica- 
tive class function associated to the polynomial P{x) = 1 — x. The main 
difference between Zn{x) and {P) is that {P) is independent of the 
interpretation of ©at as matrices. 

We now wish to obtain the asymptotic behaviour of the moments 



{W^ {Pr){x^)f^ [W\P,){x,))' 



for xi ^ X2. The easiest way to achieve this is to extend the definition of 
W^{P). 

Definition 3.4. Let P{xi,X2) be a polynomial in the two variables xi,X2. 
For a £ &N with cycle type X, we set 

e{x) 

W^{P)ixuX2) = W^{P){xi,X2)i(T) := n ^ . (3.5) 

m=l 

A simple computation using the definitions above shows 

(T^^(Pl)(xi))'^ (t^^(P2)(x2))'' = W''{P,'^){xi)W^{P^')ix2) 

= Ty^(p)(xi,x2) 

with P(X1,X2) = P^'{xi)P^%X2). 
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This shows that it is enoue; h to consider Ee [W^ {P){xi, X2)]. 

Remark. There is no restriction to the number of variables. All arguments 
used here work also for more than two variables. We restrict ourself to two 
variables since this is enough to illustrate the general case. 

3.2. Generating functions for . We compute in this section the gen- 
erating functions of the moments of multiplicative class functions. 

Theorem 3.5. Let P be a complex polynomial with 

00 

P{xi,X2)= ^ bk^^k2Xi^X2^- (3-6) 

A;i,A;2=0 



We have as formal power series 



Ee [Ty^(P)(xi,X2)] = Yl n (^^^P {x^i-,xt) (3.8) 



Y,t^hNEe[W^{P)ixi,X2)]= n {GeixI'x'^H))'"--"- (3.7) 

N=0 fci,fc2=0 

where Ge(t)^ = exp(6 • ge{t)), and Gq and ge are defined in ( p.6[ ). 
Proof. From ( |3.5D and ( |2.5| ) we have 

crebiv rn=l 

and since W'^{P) is a class function, we may use (|2.4| ) to obtain 

Ee [i^^(p)(xi,x2)] = 1- E - n ^a™^ (3.9) 

We now compute the generating function of Hj^Eq [VF^(P)(xi, X2)] with 
the help of Lemma 2.(: 

00 00 ^ ^(A) 

E t'^hr.K [W^{P){x„X2)] = E E - n ^A.^ (-1'"' 



iV=0 N=0 \\-N ^A rn=l 



X2 



A ^A 



exp E -t-'PixT, 



xf) 



m 



exp 



kl,k2=0 771=1 
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Note that 

J2 ^{x'i'x'.Hr = ge (xt^x^^t) (3.10) 

m=l 

where qq is defined in (^), and hence 



5^t^h^E[VF^(P)(xi,X2)] = n (G0(x^x^^^))''^-^2. (3.11) 

Af=0 fci,fc2=0 

and this proves (3/7), as required. □ 

Remark. The requirement for P to be a polynomial is there to ensure abso- 
lute convergence, and clearly this condition can be considerably weakened, 
see |5|, Sections 5,6]. 

As an immediate consequence we have: 
Corollary 3.6. Let si,S2 G N 6e given. Then 

j2t''hNEe[ZNixirzNix2r]= n n(G0(x^4^t))(-'^''^''(^i)(^2). 

Ar=0 ki=0k2=0 

(3.12) 

Proof. We have 

ZN{xirZM{x2y' = (iy^(l-Xi))^^(VF^(l-X2))'^ 

= w''{{i-xir{i-x2r'). 

The corollary now follows immediately by calculating the Taylor expansion 
of (1 - xiY^ (1 - X2Y^ near 0. □ 

3.3. Asymptotic behaviour of the moments. Combining the generat- 
ing functions in Theorem 3^ with the singularity analysis developed in Sec- 
tion 2^, we compute the asymptotic behaviour of Ee [M^^(P)] as N ^ oo. 

We have to distinguish between the cases \xi\ < 1 and \xi\ = 1. We consider 
here only ge{t) G T{r,'&,K). The results and computations for ge{t) G 
eJ-'{r,i),j) are similar, with only minor differences in the error terms. 

We first look at the asymptotic behaviour inside the unit disc. We have 



Theorem 3.7. Let P be as in Theorem \3.q , and let xi,X2 ^ C be given with 
max{|xi|, 1x21} < 1. Assume that g@{t) £ J^{r,-d,K), then 



with 



= E,{x,,X2) = n {Geirxl^x^^))'''-'^^. (3.14) 



(fcl,fe2)^{0,0) 
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Proof. Set 

S{t):= n {GeixI'x'^H))''-'^. (3.15) 

{fcl,fc2)^{0,0) 

Since P is polynomial, the product is finite and there is no problem with 
convergence. The domain of holomorphicity of S is thus the intersection of 
the domains of holomorphicity of each factor. This shows that the function 
S{t) is holomorphic for \t\ < r + e for an e > since max{|j;i|, \x2\} < 1 and 
GQ{t) is holomorphic for |t| < r. 

Separating the ki = k2 = term in (|3.7| ) from the rest, we can write the 
generating function as 

oo 

J^t^/iTvIEe [w''{P){xi,X2)] =exp(6o,o-5e(t))5(t). (3.16) 

N=0 

Applying Theorem |2.1(]| , we get 

..E[«.-(P)]=iV*...-'.- -l(^+o(i 

Comparing S{r) with Ei, and using Lemma [2.11| to find the asymptotic 
behaviour of /iat proves the theorem. □ 

As a special case we get the asymptotic behaviour of Ee [Z^{xi)Z^^{x2)] 
with respect to Pe inside the unit disc. 

Corollary 3.8. Let xi,X2 G C 6e given with max{|2;i|, |x2|} < 1 and let 
si,S2 G N. We then have 

Ee [Z^^{x^)Z^^^X2)] = n (Ge(rx^4^))^"'^''''''(^'^)(^'^) + O 

(3.17) 

Proof. This follows immediately from the fact that Zj\^{x) = W^{1 — x){x) 
and that (1 — xiy^{l — X2)^^ evaluated at xi = X2 = is 1. □ 

In particular, for the uniform measure {6k = 1 for all A;) Corollary 2.7 gives 
Ge{t) = (1 - t)~^ in which case we have 

E [z-(.o^^^(x.)] = n (1 - -^-r)-^^-^^^^-^^)^-^^^"^^ + o (^) . 

(3.18) 



This shows that Corollary agrees with |28, Theorem 2.13] in the unform 
case. 

The behaviour on the unit disc is more complicated. The reason is that the 
generating function can have (for fixed X\^X2) more than one singularity on 
the circle of radius r. Another point that makes this case more laborious is 
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the requirement to check whether some of the singularities of the factors on 
the RHS of ( |3.7| ) are equal. For simplicity, we assume that all singularities 
are distinct. 



Theorem 3.9. Let P be as in Theorem \3.3^ , and let xi,X2 G C be given with 
= \x2\ = 1 and x'l^X2^ / 1 for all {ki,k2) £ 1? \ {(0,0)}. Assume that 
ge{t) e J^{r,'&,K), then 

/Cl ,fv2 

with 

E2{ki,k2) = e^^^^'^i'^2-^^ n (G0(r<^^^'^x™^^^2))^'"i''"2 (3.19) 

(mi,m,2)^(fcl,A:2) 

Proof. We define 

Fit) = bk,M9e{x\'xlH). (3.20) 

fcl,fc2 

By (|3.7|) we see that exp(F(t)) is the generating function of /lArEe [VF^(P)] . 
We first take a look at the domain of holomorphicity of F{t). We have 
by assumption that gQ,{t) is holomorphic in Ao(r, i?, i;^) for a R > r and 
< (f) < ^. This shows that gQ{x^^X2^t) is holomorphic for t in the do- 
main Ai{r, R, (pjrx^''^ X2'^^) with Ai as in Definition 2.12| , and that F is 
holomorphic in 

D:= fl Ai(r,i2,</.,rx^'=ix2-'=2) = Ad(r, (3.21) 



where ^ the finite sequence of all rxJ^^X2^^ with bki^k2 ¥^ (in any order). 
Notice that this is only a finite intersection since P is a polynomial. Since 
l^^il = 12^21 = li we see that D has a shape as in figure § and that F 
has singularities at t = rx^^^X2^'^ . We thus may use Theorem p. 14 and 
therefore need to take a look at the behaviour of F near each singularity. 
We assumed that x\^x^ ^ x^^x^'^ for (777,1,777-2) 7^ {ki,k2), which implies 
that the singularities are distinct, and thus near the point rx^'^^x'^''^ , F{t) 
has the expansion 

Fit) = 6fc,,fc,791og I I I + bk,,k2K 

\l- txYx2 /r j 

+ bm„n.29e{rxT'~''x'^'~'n + 0{t-rx^'^X2''') (3.22) 
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for t —7- rxi '^^rEo ■ This shows that we can apply Theorem 2.14. Combining 



this together with Lemma [2.11 proves the theorem. □ 



Remark. For simplicity we have assumed that all the singularities are dis- 
tinct. The modification required to cope with the case when x^^Xg^ 



for some (mi,m2) ^ (^1,^2) would appear in ( 3.22 ), but techni- 
cally there is no restriction. Such a situation, with all the details written 
out explicitly, appears in p. 

To illustrate this theorem, we will calculate the autocorrelation of two char- 
acteristic polynomials at distinct points xi,X2 on the unit circle subject to 
/ foj. {fci, fca} / {0, 0}. 

The four coefficients of Zn{xi)Zj^{x2) are easy to calculate, being 60,0 = 
61^1 = 1 and 61^0 = ^0,1 = ~1 and this enables an immediate simplification 
to occur by observing that only the terms with ^fc^^^j maximal contribute; 
the others are of lower order, in this case being 0{N"'^^). 

Substituting these values into the theorem we have 
Ee [Zn{xi)Zn{x2)\ = £2(0, 0) + ^2(1, l)xf + O ( ^) +0^ ^ 



with 



N J \N^^^ 

{3.23) 

and 

E^n, 1 = ^^V^^ — ^ ^ -, (3.25 

^ Ge{rx^')Ge{rx^') ^ ' 



4. Traces 



4.1. Traces of permutation matrices. We consider in this section the 
asymptotic behaviour of traces of permutation matrices. Powers of traces 
and traces of powers have received much attention in the random matrix 
literature (see for instance Q, More specifically, we first look at 

Tr(o"'^) for fixed d G Z. Since the embedding of ©at into the unitary group 
in ( |3.lD is a group homomorphism, we can interpret as d-fold matrix 
multiplication and as the matrix corresponding to g o • — o . 

d times 

We first recall a well known explicit expression for Tr(cj'^) that we shortly 
prove for completeness. 

Lemma 4.1. We have for d £ Z 

TV 

Tr{a^) = Y,Mdf'Ck{a), (4.1) 
fc=i 
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Proof. The matrix corresponding to cr'^ has the form o-d(j)). We thus have 

TV 



I¥(a'^)=^<5,,,.(,)=#{i:a'^(i)=i}. 



(4.2) 



8=1 



Therefore Tr(c7'^) is the number of 1— cycles of a'^. A simple computation now 
shows that the number of 1— cycles of a'^ is indeed Ylk=i lfc|d^Cfc(o')- ^ 

Using this expression and the method of generating functions developed in 
Section 2^, we prove a weak convergence result for Tr(cj'^). 

Theorem 4.2. Let d £N be given. We then have 



J2 



e 



N=0 



t^ = e^vij2jie'''-ir]Geit). 

k\d 



(4.3) 



If 9e is of class J^{'d,r, K), then 







Js Tr(o-'* 



exp 



k / Jsk 



k\d 



l)r^ +0 



(4.4) 



If 9e is of class eT{'0,r,^), then 



E© 



exp 



^^(e-'=-l)rM +0 



k\d 



log(iV) 

N'T 



(4.5) 



Proof. Applying Lemma [4.1| , and evaluating the expectation explicitly in 
terms of partitions using Lemma ^.3| , we have 

oo 



e 



N=0 



X ^ k=l 



Cfe 



(4.6) 



The cycle index theorem (Lemma 2.6) yields that this equals 



exp 



/ oo ^ 



\k=l 



^^gi«fci,|,^fc =exp 



ET(^"'-l)*'|Ge(t) (4.7) 

k\d 



where G0(t) is given in (|2.6| ). This proves equation ( [4.31) . 

Applying Theorem 2.10| to this yields equation ( |4.4| ), and Theorem 2.16 
yields equation (|4.5[), as required. □ 



Remark. An alternative way to prove Theorem is to use Theorem 3.1 in 
[23|, which computes the generating function of /lArEe exp (Yl^k=i i^kCk^ 
and its asymptotic behaviour for ge{t) G F{'&,r,K) and ge{t) G eJ^('!?, r, 7). 

We obtain as an immediate corollary: 
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Corollary 4.3. Let d £ 7, be fixed and assume that qq is in r, K) or 

eJ^(??, r, 7). Then 

Tr(cj'^) -^^kPk as N ^00 (4.8) 

k\d 

where are independent Poisson distributed random variables with E [P^] = 
k ^ ■ 

4.2. Traces of functions. Recall from the introduction that if M is the 
permutation matrix representing the permutation a, then for a function 
F : S*^ — )• C , we defined the trace of F to be the function Tr(F) : ©at — )• C 
with 

N 

Tr{F){a):=^F{u;k) (4.9) 

k=l 

where {ujk)k=i eigenvalues of M or o" with multiplicity. Moreover, 

N 

Tr{F){a) = ^kCkAk{F) (4.10) 

k=l 

with 

1 

Afe(F) :=-^F(e2^^'"/'=) (4.II) 

m=l 

(compare with Lemma |l.2| ) . 

The asymptotic behaviour of Ti{F) is not so easy to compute for an arbitrary 
function defined on the unit circle. This problem will be dealt with more 
carefully in Sections ^ and ^. However, if F is a Laurent polynomial, we can 
use the same method as for Tv{a'^). 

Theorem 4.4. Let 

F{x) = Y,bdx'' (4.12) 

d 

be a Laurent polynomial. If ge € J-{'0,r,K) or gQ G eJ-'(??, r, 7), then 

00 

Tr(F)(cj) -Nbo^ bdY^ kPk as N ^ 00 (4.13) 

d=—oo k>l 
d^O k\d 

where Pj. are independent Poisson distributed random variables with E [P^] = 



22 



C. HUGHES, J. NAJNUDEL, A. NIKEGHBALI, AND D. ZEINDLER 



Proof. Due to the linearity of Tr(F), we may assume the constant term, bo, 
is zero. As before we apply the cycle index theorem to obtain 







N=0 



JsTr(F)ia) 



A k=l 



exp 



E 

\k=l 



(4.14) 
(4.15) 



and since A^^x'^) = and is linear, this equals 



,k=l 



exp j ^ ^ exp I is ^ bdk l^aj 



exp 



Et n 



A;=l 



v 



Jsb^k 



\ k\d 



I 



Geit) (4.16) 



Note the first factor is entire, so Theorem [2.10| (for the case of J-{r,'d, K)) 
and Theorem 2.16| (for the case of eJ-{r,'d,j)) yields 



E. 



e 



exp 



E 

k=l 



n 

\ k\d 



^isb^k _ 2 



k 



(4.17) 



The RHS is the characteristic function of the RHS in ( [4.13| ) . This completes 
the proof. □ 

Theorem 4.5. Let F : ^ C be of bounded variation, and d G N 6e given. 
Then 



1 



Tr(F)(a; 



I^J{v)d^^ +0 



N 



(4.18) 



where T{a) is the total number of cycles of a, and dip the uniform measure 
on . 

Moreover if geit) G J^{i},r,K) or gQ{t) £ eJ^(i9,r, 7) then Ee [T(cr)] ~ 
'!91og(A^), and thus we have a quick convergence of the moments. 



Proof. Since F is of bounded variation we can apply Koksma's inequality 



(|19, Theorem 5.1]) to see that 



k 

m=l •^•^^ 



dip 



< 2DkV{F) 



(4.19) 
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with V{F) the variation of F and the discrepancy of the sequence 
^g27rim/fc'j^^^^ But the discrepancy is dominated by 1/k. We thus have 

^^(^) = ^ E F{e'-'^/') = I F{^) + (i) (4.20) 

m=l ''^^ ^ ^ 

We now combine ( |4.1(]| ) and ( [4.20| ) and get 

N 



Tr{F)ia) = j2Ck(k[ Fi^p) dip + 0{r 



= N / F{ip) d^ + 0{T{a)), 

(4.21) 

where we have used that X^^Li kCk = N and '^^=1 = T{(j). Notice that 
( |4.21| ) is independent of any probabihty measure on ©at. Using the binomial 
theorem and the fact that < T{a)/N < 1 for all a, we get 

^ {TtiF){a)r = (1^^ F(^) dy.y + O,,, (^) (4.22) 

where the constant implicit in the big-0 is independent of a and N. We 
apply Ee [.] on both sides, and this proves the first part of the theorem. 

The last statement follows from Theorem 4.2], where it is shown that if 
geit) G F{^,r,K) or ge{t) G eJ-(i9,r,7) then Eq [T(a)] ~ i}log{N). □ 

Remark. In fact, for many probability distributions on &n, Ee [r(c7)] = 
o{N). The only way for this not to be true is for a to frequently have only 
small cycles, which will occur if = {0k)^=i is a sequence tending to zero 
very rapidly. 

5. Wreath product, traces and the generalised Ewens measure 

We consider in this section the traces of the wreath product 16n (see for 
instance ||2^ ) . More precisely, we consider random matrices of the form 

M{a, zi, zjy) := diag(2;i, . . . , zat) • a 

where a is a random permutation of &n, and {zj)j>i is a sequence of i.i.d. 
random variables with values in (the complex unit circle), independent 
of a. Many groups closely related to ©tv give similar matrices, for instance 
the Weyl group of SO (27V). 

The trace of a function F is then extended in the obvious way by 

N 

Tr(F) =TV(F,z)(a) ■.= Y,H^k) (5.1) 

k=l 

where {(^k)k=i ^ eigenvalues of M{a,z). 

We now give a more explicit expression of Tr(i<'). 



24 C. HUGHES, J. NAJNUDEL, A. NIKEGHBALI, AND D. ZEINDLER 

Lemma 5.1. 

N Ck 

Tr(F) = Y,Y. ^^k{F, Zk,m) 

k=l m=l 



(5.2) 



where {Zk^m)k,m>i is a sequence of independent random variables which is 
independent of {Ck)k>i (the sequence of cycle numbers of a), with Z^^m equal 
in distribution to 11^=1 > '^^'^ 



(5.3) 



Proof. The characteristic polynomial of M{a, zi, zn) with a £ &N with 
cycle type A, is given by 



N Ck 



det (1 - xM{a, zi, zn)) = H 11 U " 11 4 



k,m 



(5.4) 



k=l m=l 



where the sequence {z^j'^)k,m,j is the same sequence as (-Zj)jLi) but with a 
different numeration and ordering. (Note that this is why it is crucial that 
the {zj) are i.i.d.). The proof of ( |5.4| ) is similar to the proof of ( |3.3|) a nd we 
thus omit the details. The lemma now follows immediately from ( |5.^ . □ 

As in the last section, we can compute the generating function of Tr(F). 
Lemma 5.2. We define 



Xk{s) :=]E 



JskAkiF,Zk,m) 



We then have 







JsTr{F) 



(5.5) 



(5.6) 



N=0 



\k=l 



Remark. Note that Xfc(s) is independent of m since Z^^i = Zk^2 = 



Z. 



k,m • 



Proof. We compute E [exp(zs Tr(F))] . For this we use the independence of 
Ck and Afc to obtain 



E 



,is Tr{F) 



E 



■ N Ck 
,fc=l m=l 



JskAkiF,Zk,m) 



E 



E 



■ N Ck 

n n xkis] 

,k=l m=l 

■ N 

.fc=i 



The theorem now follows immediately from Lemma 2.6 



(5.7) 

□ 
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Definition 5.3. Let 



oo 



k 

k=l 



Theorem 5.4. Assume E[\Ak{F, Zk,i)\] = 0(A;-^-^) for some < 5 < 1, 
and assume ge{t) is in eJ-"(r, t?, 7), where eJ-"(r, t?, 7) is given in Defini- 
tion 2.1 1 . Then 5'Tr{F)(*) eJ^(r, min {7, (5}) and 



gisTr(F) 



«^p(Et(X'^(^)-1)^') +0(iV-"'''^^^'^hog(iV)) (5.9) 



\k=l 



and as tends to infinity, Th:{F) converges in law to the random variable 

00 Pfc 

Y:=Y,Y. ^^k{F, Zfc,„), (5.10) 

k=l m=l 

where {Pk)k>i is a sequence of independent Poisson random variables, inde- 
pendent of {Zk,m)k,m>i, o,nd such that has parameter O^r^ /k. Here, the 
series defining Y is a.s. absolutely convergent. 

Remark. By linearity of trace, if F is Riemann integrable one can always 
subtract a suitable constant to make F{(p) d(p = 0, which ensures 
Afc(F, z) — )• as A; — )• 00. 

Remark. One should compare equation ( 5.1C| ) with equation ( |5.2| ). The 
replacement of the cycle counts with P^ is indicative of Feller-coupling 
for the generalised Ewens measure. 



Proof. We have 

00 „ 00 „ 

9TriF){t) = -fxk{s)t'^ = ge{t) + ^ -±{xk{s) - l)t'= (5.11) 
k=l ■ ' 



k=l 



We now have 

\Xkis)-l\ <E 



JskA^,{F,Zk,l) _ I 



<E[ik\AkiF,Zk,i)\)] = Oik- 



(5.12) 



On the other hand we have Ok = 0{r ^). This follows immediately from the 
fact that ge{t) is in eJ^(r, i?,7). We thus have 



{Xk{s) - 1) = 0{r-^k-^-'). 



(5.13) 



This shows that fi'Tr(i^)(*) ^ eJ^(r, i9, min {7, (^}). 
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Since ge{t) G eJ"(r, ??, 7), we can write ge{t) = i?log(^Y^^j + go{t) with 
go{r) < 00. Thus 

5Tr(F)(t) = ^log (y^) +9o{t) + p^^f{Xk{s) - l)tK (5.14) 



We get with Theorem [2.16| that 

Js Tr{F) 







N'' 



exp 50 (r) + -^(Xfe(s) - + Rn 



k=l 



with 



Rn = 



(5.15) 



(5.16) 



Dividing by /iat proves equation (^). 

Using the characteristic function of Tr(F), we can deduce its convergence in 
law to Y . The absolute convergence of the series in ( ^.lOD comes from 

00 Pk 



E 



.fc=l m=l 



J^fcE [Pfc]E[|Afc(F,Zfe,i) 

k=l 



(5.17) 



k>l 
< 00, 



since 0^ = 0{r Now, for s € M and A; > 1, 



E 



>fcE™iiAfe(F,Zfe,„) 



E 



E 



JskAk(F,Zk,i) 



E 



(X.(«))'''= 



(5.18) 



ew{"-^iXk{s)-iy], 



Hence by absolute convergence, 



E [, 



AsYl 



exp 



J2iixkis)-iy 



(5.19) 



and thus by equation (|5]|), Ee [e*'^'(^)] E [e^*'*'] as iV ^ 00, and thus 
Tr(F) converges in law to Y. □ 



With a more direct approach one can prove the convergence in law of Tr(F) 
to Y (albeit without a rate of convergence) under slightly weaker conditions. 
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Theorem 5.5. Assume that ge{t) is in eJ^{r,i^,-f) and that 



k^'^'^^^E [\Ak{F, Zk,i)\] < oo. (5.20) 

k=l 

Then Tt{F) converges in law to Y, where Y is given by ( 5.10| ). 



Proof. Under these conditions, the absolute convergence of the series defin- 
ing Y is checked as follows: 



E 



oo I^k 



,k=l m=l 



J2kE[Pk]E[\Ak{F,Zk,i)\] (5.21) 

k=l 

oo 

Y,Okr''E[\Ak{F,Zk,i)\] (5.22) 

k=l 

o(f;E[|Afc(F,Zfc,i)|] ) (5.23) 



\k=l 



and this converges by assumption. 

In ||2^, Corollary 3.1.1] it is proven that for all fixed b > I, (Ci, C2, . . . , Cj,) 
tends in law to {Pi,P2, ■ ■ ■ , Pb) when the dimension goes to infinity. 

Now let 



and 



b c\ 

Trb(F) :=^^ fcAfc(F,Zfc, 

fc=l m=l 



b Pk 

n:= J] J^fcAfc(F,Zfc,^) 

k=l m=l 



(5.24) 



(5.25) 



The same argumentation as in Theorem 5.4 gives 



t^ = exp ( ^ ^fiXkis) - l)t'' ) e^eW (5.26) 



Af=0 \fc=l 

from which follows (again by the same reasoning as in Theorem 5.4) 



E [e 



(5.27) 
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as — )• oo. Here b is fixed but arbitrary, so this convergence implies (by 
using the inequality |e*^ — e*^| < |x — ?/|) that 



hm sup 

7V-!>oo 



JsTt{F) 



E [, 



isYl 



<|s|hmsup J2 kK[\AkiF,Zk,i)\]E'[iCk + Pk)] 
k=b+i 

oo 

< |s| J2 m\^k{F,Zk^^)\]Hk (5.28) 



k=b+l 

where E' is the expectation over the product measure of Pe and the measures 
occurring from {Pk)k>i, and where 

= E[Pk] + sup Ee [Ck] . (5.29) 

7V>1 

Therefore the theorem is proven if we show that 

oo 

Y^kHkE [|A,.(F,Zfc,i)|] <cx). 

k=l 

Ercolani and Ueltschi Proposition 2.1(c)] show that 
Ee [Ck] 



(5.30) 



k /ijv 

if A; > 



(5.31) 



By Lemma p. 17 , we have, for some A> and for N going to infinity, 

hN^A{N + lf~^/r^, (5.32) 

so 

hN~k 



Now, for k fixed. 



max I 1 

N>k 



k 



hN 



O rM 1 



k 

N + 1 



(5.33) 



1 if 1? > 1 

{k + 1)^-^ if 1? < 1 (attained at iV = A:). 

(5.34) 



Since g@{t) G eF{r,-d,'^) we have O^r^ = 0(1), and so we deduce that 

sup Ee[Cfc] = of^r'= (A: + 1) ^ = 0{k-^+^-^-'^^+). (5.35) 

Ar>l V ^ / 

Finahy, since 



we have 



E [Pk] = Okr'/k = 0{l/k) 
H, = 0{k-'+^'-^H 



(5.36) 
(5.37) 
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oo / oo \ 

^kHkE[\Ak{F,Zk,i)\]=oi^k^'~^^+K[\Ak{F,Zk,i)\]\ < oo (5.38) 

k=l \k=l J 

as required. □ 

Remark. In the case when Zi are all equal to 1 almost surely, then Afc(F, = 
Ak{F) as given in ( [4.11| ), and we are back in the case of permutation ma- 
trices. Thus these two theorems fulfill the promise made in Section ^. 

The following result gives sufficient conditions, expressed only in terms of 
the Fourier coefficients of F, under which we can be assured the conditions 
of Theorems 5^ and |5]^ apply. 

Theorem 5.6. Let us suppose that F is continuous and for m G let us 
define the m-th Fourier coefficient of F by 

c^{F) := — /'^^e-*™^F(e^^)(ix. (5.39) 
271" Jo 

We assume that the mean value of F vanishes, i.e. cq{F) = 0. If for some 
5 € (0, 1], Cm{F) = 0{\m\^^~^) when \m\ goes to infinity then 

E[\AkiF,Zk,i)\] = 0{k~'-^). (5.40) 

If there exists s > (1 —"!?)+ such that 

\mncmiF)\ < oo (5.41) 

meZ 

then 

oo 

k^'-^^-^E [\Ak{F, Zk,i)\] < oo. (5.42) 

k=l 

Remark. If the assumptions of Theorem are satisfied, except that cq{F) = 
0, then one can still apply the result to the function F — co{F), and deduce, 
from Theorem or Theorem 5.5, that Tr(i^) — Nco{F) converges in law 



to Y, where Y is given by (5.10). 



Proof. Since F is continuous, one has, for all x G [0, 2tt), 

Fie''') = lim y i!i^-M±c^(i7)e*™^, (5.43) 

n— >oo ^—^ n 



by using the Fejer kernel. Now, by assumption, 

Y \Crn{F)\ < OO, (5.44) 

and hence, by dominated convergence, 

F(e*-) = Y cm(i^)e'"^^ (5.45) 

meZ 
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where the series is absolutely convergent. Since co{F) = 0, one deduces that 
for all > 1 and x £ [0, 27r), 



j=0 



k-l 



meZ\{0} \j=o 

meZ\{0}, 
k\m 

If F satisfies the first assumption, that Cm{F) = 0(|m|"^"^), then 

/ \ 



sup 

xe[0,27r) 



|Afc(F,e^^)|< Yl \Cra{F)\=0 



meZ\{0}, 
k\m 



E 

mGA{0}, 
k\m 



\m\ 



(5.46) 



(5.47) 



for k going to infinity, which clearly implies E [| Afc(F, = 0{k ^ ^) 

If F satisfies the second assumption, one has 



sup |Afc(F,e-)l <5]A;(i-'')+ ^ |c„(F)| (5.48) 

k\m 



k=l 



<E E \^\^^~''^^\cm{F)\ (5.49) 



k=l meZ\{0}: 
k\m 

< J2 |c™(i^)lkl^'"''^+T(|m|), (5.50) 

meZ\{0} 

where T(|m|) denotes the number of divisors of \m\. Since T(|m|) = Odml*^) 
for all e > 0, one deduces that 

J]A;(i-'')+E[|Afe(F,Zfc,i)|] = I ^ \cm{F)\\mA < ^. (5.51) 
fc>l \m.ez\{o} / 

This completes the proof of the theorem. □ 
One deduces the following: 

Corollary 5.7. Let F be a continuous function from to C, contained in 
a Sobolev space for some s > 1/2 + (1 — "!?)+. Then, the second condition 



of Theorem 5.t is fulfilled, and thus also the conditions of Theorem 5.5 
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Proof. By Cauchy-Schwarz inequality, one has, for any a E ((1 — '&)+,s — 
1/2), 

1/2 / \ 1/2 



meZ\{0} \meZ\{0} / \meZ\{0} 

(5.52) 

which is finite since F G and 2(q — s) < —1. □ 

Remark. Note that it is not always obvious to estimate directly the Fourier 
coefficients of a function F, however standard results from Fourier analysis 
concerning the differentiability of F yield sufficient bounds on the decay of 



the Fourier coefficients of F for the conditions of Theorem 5^ to be checked. 
See, for example, [|^, Chapter 9]. 

6. Diverging variance for the classical Ewens measure 

In the previous two sections, we've been considering the convergence of 
Tr(F) to some limit for random permutation matrices (and their general- 
isation to wreath products), where the underlying probability space is the 
generalised Ewens measure. The conditions we've used have all implied that 
the variance of Tr(i<') stays bounded as N ^ oo. 

A recent paper by Ben Arous and Dang ||^] dealing with Tv{F) for random 
permutation matrices in the special case of the classical Ewens measure 
(which is when 9k = 0, a. constant), demonstrates a dichotomy between con- 
verging and diverging variance for Tr(F) in the classical Ewens distribution. 
In the former case they also show convergence of Tr(F) to an explicit finite 
limit, and in the latter case they prove the following central limit theorem: 

Theorem 6.1 (Ben Arous and Dang). Assume that 

N 

Vn := eY,kAk[Ff (6.1) 

k=l 

tends to infinity as N ^ oo and 



then 



^max^fc|Afc(F)| =o( V^^) (6-2) 



Tr(F) - E [Tr(F)] d , , , 

^ ^ L ^ ^ ^ AA(0,1). (6.3) 



VVn 

In the generalised Ewens measure, we are currently unable to apply the 
function theoretic methods to prove weak convergence results in the case of 
diverging variance. However, for the classical Ewens measure we are able 
to prove a similar central limit theorem the wreath product, with slightly 
extended application, in the sense that condition ( |6.2| ) can be weakened from 
a sup-norm to a p-norm. 
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Theorem 6.2. Assume that 

N 

Vn ■.= eY,k^ [{^k{F,Zk,i)f\ (6.4) 

k=l 

tends to infinity as N — )• oo. Assume further that there exists a p > 
max {^,2 1 such that 

N 



^kP-'E[\Ak{F,Zk,i)n = o(v^/') (6.5) 

k=l 

with Ak{F, z) = i Ylu^.^z ^(^)- Then 

I^Tr(F) - En 



J ATM 



converges in distribution to a standard gaussian random variable, where 



N 



EN:=eyE[Ak{F,Zk,i)]. (6.7) 



k=l 



Remark. Recall that without loss of generality we may assume F has mean 
zero in the sense that J^^ F{e^^)dx = 0. We remark that this does not 
necessarily imply that Afc(F, z) tends to zero, even though Afc(F, z) is a 
discretisation of the integral, without the assumption of further smoothness 
conditions. Moreover, note that in the framework of the symmetric group 
(i.e. Zfs i = 1 almost surely), the assumption ( |6^ ) is implied by the condition 
( |6.2[) given in Indeed, if ( |6.2[) is satisfied, then for p > 2, one has 

^F"i|A,(F)r < f max^(A;|A,(F)|) j;fc|Afc(F)|2 

k=l ^ ^ k=l 

= o{vt'^/')0{VN) = o{vf). 

In the proof of Theorem |6.2| , we use Feller coupling, which allows the random 
variables Ck and Pk to be defined on the same space. This coupling allows 
us to prove the theorem with replaced by P^, and the following lemma 
allows us to estimate the distance between the two. 

Lemma 6.3 (Ben Arous and Dang [||). For any 6 > there exists a con- 
stant K[9) depending on 9, such that for every 1 < m < N, 

E[\Ck-Pk\]<^ + ^^N{k), (6.8) 



where 



'N - k + e -1\ fN + 9 - ^ , , , 
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Proof of Theorem 6A . The main idea of the proof is to define the auxihary 
random variable 



N Pk 

Yn{F) ■.= ^^kAkiF,Zk,m) 



.10) 



k=l m=l 



and to show that Tr{F) and Yi'^{F) have the same asymptotic behaviour 
after normahsation, and that (again after normahsation) Y/v(F) satisfies a 
central limit theorem. 

First we will show that 

E[\Ti{F)-Yn{F)\] = o(^{Vn)'^) (6.11) 



We use Lemma 16. 31 and that Z/^ 

are independent of Ck and Pk to get 



E[|Tr(F)-y^(F)|] <E 



N 



Y,\Ck-Pkm[\Ak{F,Z^ 



k,l) 



.k=l 



< ^ E ^'^m)!] + nY.^^ ^'^m)!] ^^(^) (6.12) 

k=l k=l 



For the first term, we apply Jensen's inequality and condition (|6.5|) to obtain 

N / N \Vp 

-Y.m[\Ak{F,Zk,i)\]< -J]FE[|Afe(F,Z,,i)|f (6.13) 



k=l 



k=l 



< (^5^F-iE[|A,(F,Z,,i)n^ 



i/p 



o V, 



1/2 



TV 



(6.14) 
(6.15) 



Now we deal with the second term in ( 6.12 ). If ^ > 1 then ^'jv(/c) is bounded 
by 1, so the same argument as above shows that the second summand is also 



in this case. 



If < < 1, we have to be more careful. A simple computation shows that 
there exists constants Ki,K2 such that 



y^7v(fej<i j^^^i-e^ for k = N 



.16) 
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and so 
^ 

-Y,kn\^k{F,Zk,imN{k) (6.17) 



N 

k=l 



k=l ^ 



Using the value of p given in the conditions of the theorem, 



(6.18) 



N'-'E[\An{F,Zm,i)\] = [NP-P'E[\AN{F,Z^^i)\ry^'' (6.19) 

< ( j;F-P'E[|A,(F,Zfc,i)|H (6.20) 

< ( J^F-iE[|Afc(F,Zfc,i)r]j (6.21) 



since kP'P'' < k^'^ (since pO > 1) and E [\Ak{F, Zk,i)\]^ < E [\Ak{F, Zk,i)\'P] 
(since p > 1). Thus, by condition (6^) this is o (^V^^ 

Holder's inequality gives 



^^kE[\AkiF,Zk,,)\](l 



k=l 



<[-X:FE[|A.(F,Z,,)nj (^^E(l-]^) ) (6-22) 



with 1 + 1 = 1. 

V q 



After re-ordering the sum, the second factor is [a^'^*-"'^ ^Yl!j=i j '^^^ ^^'^ 



■'J 

and if q{6 — 1) > — 1 then it is bounded above by a constant. Note that 

{e-i)q>-i ^ {i-e)<- ^ (1-0) <i-- ^ e>- 

q p p 

and thus condition ( |6.5| ) now ensures the existence of a p > ^ such that 
the first factor is o ^(y/v)^^ and the second factor is bounded. This proves 

(iH). 

Therefore Slutsky's theorem implies that Yn{F) and Tr(i<') have the same 
asymptotic distribution after scaling. Thus it suffices to show that 

ME^^ (0.23) 
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converges in law to a standard gaussian random variable. 

We calculate the mean of Y/v(F) by first taking expectation with respect to 
Zfc^m and then with respect to Pk, to obtain 



N 

= ^E[Pfc] kE [Ak{F,Zk,i)] 

k=l 

(6.24) 

where we use the fact that E [Afc(F, Zf^^m)] = E [Af^^F, .^fc,i)] for all m. Fi- 
nally, since E [Pk] = 6/k, we see that E [Yn{F)] = En as defined in (|6J|). 

For the variance, since Pk and Zk^m are all independent, one can move the 
sum outside the variance, to obtain 



N 



E[y^(F)] = J^E 



k=l 



Pk 

E 

_m=l 



kE [Afc(F, Zk,m)] 



N 



Var(y,v(i^)) = J^Var kAk{F,Zk,r 



3.25) 



k=i 



\m=l 



Now, the variance of a sum of random length of i.i.d. random variables is 
given by the following formula: 



Var ( X] ^™ 1 = Var(Xi)E [P] + Var(P)E [Xi 



(6.26) 



\m=l 



if {Xm)m>i are i.i.d., random variables, and if P is a variable, in- 
dependent of {Xm)m>i (this result can be proved by a straightforward cal- 
culation). Letting X^ = kAk{F, Zk^m) and P = Pk and knowing that 
E [Pk] = Var(Pfc) = 0/k, we deduce that Var(yAr(F)) = Vn- 

Finally we apply the Lyapunov central limit theorem since Yn{F) is a sum 
of independent random variables. We will show that 



N 



k=l 



Pk 



m=l 



kAk{F, Zk,m) - E [kPkAk{F, Zk,i] 

N 

«5^F-iE[|Afc(F,Zfc,i)|^] (6.27) 



fe=i 



and by Condition ( |6.5D with p > 2 this is o{V^ ) which means ^'^^y 
converges in law to a standard gaussian random variable. 

For simplicity, let P be a Poisson random variable with parameter 9/k (we 
think of k as being large), and let Xm = kAk{F, Zk^i) be i.i.d. random 
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variables with E finite. Then 



E 



m=l 



X„ - E [P] E [Xi] 
p 

Y,{Xm-^ [Xl]) 



< E 



m=l 



E 



Y,iXm-^ [Xl]) + (P - E [P])E [Xl] 



m=l 



+ E[\P-E[P]\PfPE[Xi]\ (6.28) 



by the generahsed triangle inequality. 
Now, for all p > 1, 

E [|P - E [P] \P] < E [PP] < E [P] 
as E [P] — 5- 0, and so the second term in ( |6.28D is 

E [|P - E [P]|P]^/pE [Xl] < E [P]^/^'E [Xl] 

< (E [P]E[|Xi|P])'/^ 

by Holder's inequality, since p > 1. 



(6.29) 

(6.30) 
(6.31) 



To bound the first term in ( |6.28D , let g„ = P [P = n] , and note that 



E 



m=l 



n=0 

oo 



5^(X^-E[Xi]) 



m=l 
n 



<^(ZJ 5;E[|X^-E[Xi]n 



pii/p 



71=0 \m=l 

oo 



^g„nfE[|Xi-E[Xi]|P] 



n=0 



= E[pP]E[|Xi -E[Xi]|P] 
<E[P]E[|Xi|P] 



(6.32) 
(6.33) 

(6.34) 

(6.35) 
(6.36) 



Thus 



E 



^X^-E[P]E[Xi] 



m=l 



« ( (E [P] E [|Xi|P])^/^' + (E [P] E [\Xi\P])'^Py 



<E[P]E[|Xi|P] 



(6.37) 



Using E[P] = e/k and E[|Xi|p] = kPE[\Ak{F, Zk,i)\^], and summing for 
A; = 1, . . . , iV we have proven ( |6.27| ). By Condition (|6.5|), if p > 2 then this 
is o(V^^^) which by Lyapunov's theorem means that XEi£^MlL converges in 
law to a standard gaussian random variable as required. □ 
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